Constraining the cometary flux through the asteroid belt during the late heavy bombardment

In the Nice model, the late heavy bombardment (LHB) is related to an orbital instability of giant planets which causes a fast dynamical dispersion of a transneptunian cometary disk. We study effects produced by these hypothetical cometary projectiles on main-belt asteroids. In particular, we want to check whether the observed collisional families provide a lower or an upper limit for the cometary flux during the LHB. We present an updated list of observed asteroid families as identified in the space of synthetic proper elements by the hierarchical clustering method, colour data, albedo data and dynamical considerations and we estimate their physical parameters. We selected 12 families which may be related to the LHB according to their dynamical ages. We then used collisional models and N-body orbital simulations to gain insight into the long-term dynamical evolution of synthetic LHB families over 4 Gyr. We account for the mutual collisions, the physical disruptions of comets, the Yarkovsky/YORP drift, chaotic diffusion, or possible perturbations by the giant-planet migration. Assuming a "standard" size-frequency distribution of primordial comets, we predict the number of families with parent-body sizes D_PB >= 200 km which seems consistent with observations. However, more than 100 asteroid families with D_PB >= 100 km should be created at the same time which are not observed. This discrepancy can be nevertheless explained by the following processes: i) asteroid families are efficiently destroyed by comminution (via collisional cascade), ii) disruptions of comets below some critical perihelion distance (q <~ 1.5 AU) are common. Given the freedom in the cometary-disruption law, we cannot provide stringent limits on the cometary flux, but we can conclude that the observed distribution of asteroid families does not contradict with a cometary LHB.Comment: accepted in Astronomy and Astrophysic

Did the Hilda collisional family form during the late heavy bombardment?

We model the long-term evolution of the Hilda collisional family located in the 3/2 mean-motion resonance with Jupiter. Its eccentricity distribution evolves mostly due to the Yarkovsky/YORP effect and assuming that: (i) impact disruption was isotropic, and (ii) albedo distribution of small asteroids is the same as for large ones, we can estimate the age of the Hilda family to be $4_{-1}^{+0}\,{\rm Gyr}$. We also calculate collisional activity in the J3/2 region. Our results indicate that current collisional rates are very low for a 200\,km parent body such that the number of expected events over Gyrs is much smaller than one. The large age and the low probability of the collisional disruption lead us to the conclusion that the Hilda family might have been created during the Late Heavy Bombardment when the collisions were much more frequent. The Hilda family may thus serve as a test of orbital behavior of planets during the LHB. We tested the influence of the giant-planet migration on the distribution of the family members. The scenarios that are consistent with the observed Hilda family are those with fast migration time scales $\simeq 0.3\,{\rm Myr}$ to $3\,{\rm Myr}$, because longer time scales produce a family that is depleted and too much spread in eccentricity. Moreover, there is an indication that Jupiter and Saturn were no longer in a compact configuration (with period ratio $P_{\rm S}/P_{\rm J} > 2.09$) at the time when the Hilda family was created

Shape models and physical properties of asteroids

Despite the large amount of high quality data generated in recent space encounters with asteroids, the majority of our knowledge about these objects comes from ground based observations. Asteroids travelling in orbits that are potentially hazardous for the Earth form an especially interesting group to be studied. In order to predict their orbital evolution, it is necessary to investigate their physical properties. This paper briefly describes the data requirements and different techniques used to solve the lightcurve inversion problem. Although photometry is the most abundant type of observational data, models of asteroids can be obtained using various data types and techniques. We describe the potential of radar imaging and stellar occultation timings to be combined with disk-integrated photometry in order to reveal information about physical properties of asteroids.Comment: From Assessment and Mitigation of Asteroid Impact Hazards boo

Constraining the Physical Properties of Near-Earth Object 2009 BD

We report on Spitzer Space Telescope IRAC observations of near-Earth object (NEO) 2009 BD that were carried out in support of the NASA Asteroid Robotic Retrieval Mission (ARRM) concept. We did not detect 2009 BD in 25 hrs of integration at 4.5 micron. Based on an upper-limit flux density determination from our data, we present a probabilistic derivation of the physical properties of this object. The analysis is based on the combination of a thermophysical model with an orbital model accounting for the non-gravitational forces acting upon the body. We find two physically possible solutions. The first solution shows 2009 BD as a 2.9+/-0.3 m diameter rocky body (rho = 2.9+/-0.5 g cm-3) with an extremely high albedo of 0.85(+0.20/-0.10) that is covered with regolith-like material, causing it to exhibit a low thermal inertia (Gamma = 30(+20/-10) SI units). The second solution suggests 2009 BD to be a 4+/-1 m diameter asteroid with pV = 0.45(+0.35/-0.15) that consists of a collection of individual bare rock slabs (Gamma = 2000+/-1000 SI units, rho = 1.7(+0.7/-0.4) g cm-3). We are unable to rule out either solution based on physical reasoning. 2009 BD is the smallest asteroid for which physical properties have been constrained, in this case using an indirect method and based on a detection limit, providing unique information on the physical properties of objects in the size range smaller than 10 m.Comment: 28 pages, 8 figures, accepted for publication in Ap

Stress field and spin axis relaxation for inelastic triaxial ellipsoids

A compact formula for the stress tensor inside a self-gravitating, triaxial ellipsoid in an arbitrary rotation state is given. It contains no singularity in the incompressible medium limit. The stress tensor and the quality factor model are used to derive a solution for the energy dissipation resulting in the damping (short axis mode) or excitation (long axis) of wobbling. In the limit of an ellipsoid of revolution, we compare our solution with earlier ones and show that, with appropriate corrections, the differences in damping times estimates are much smaller than it has been claimed. This version implements corrections of misprints found in the MNRAS published text.Comment: 14 pages, 6 figures, published in Monthly Notices RAS (containing misprints

New inclination changing eclipsing binaries in the Magellanic Clouds

Context: Multiple stellar systems are unique laboratories for astrophysics. Analysis of their orbital dynamics may reveal invaluable information about the physical properties of the participating stars. Unfortunately, there are only a few known and well described multiple systems, this is even more so for systems located outside the Milky Way galaxy. A particularly interesting situation occurs when the inner binary in a compact triple system is eclipsing. This is because the stellar interaction, typically resulting in precession of orbital planes, may be observable as a variation of depth of the eclipses on a long timescale. Aims: We aim to present a novel method to determine compact triples using publicly available photometric data from large surveys. Here we apply it to eclipsing binaries (EBs) in Magellanic Clouds from OGLE III database. Methods: We analyzed light curves (LCs) of 26121 LMC and 6138 SMC EBs with the goal to identify those for which the orbital inclination varies in time. Archival LCs of the selected systems, when complemented by our own observations with Danish 1.54m telescope, were thoroughly analyzed using the PHOEBE program. Time dependence of the EB's inclination was described using the theory of orbital-plane precession. By observing the parameter-dependence of the precession rate, we were able to constrain the third companion mass and its orbital period around EB. Results: We identified 58 candidates of new compact triples in Magellanic Clouds. This is the largest published sample of such systems so far. Eight of them were analyzed thoroughly and physical parameters of inner binary were determined together with an estimation of basic characteristics of the third star. These data may provide important clues about stellar formation mechanisms for objects with different metalicity than found in our galactic neighborhood.Comment: Accepted for publication in Astronomy and Astrophysic

Candidates for asteroid dust trails

The contribution of different sources to the circumsolar dust cloud (known as the zodiacal cloud) can be deduced from diagnostic observations. We used the Spitzer Space Telescope to observe the diffuse thermal emission of the zodiacal cloud near the ecliptic. Several structures were identified in these observations, including previously known asteroid dust bands, which are thought to have been produced by recent asteroid collisions, and cometary trails. Interestingly, two of the detected dust trails, denoted t1 and t2 here, cannot be linked to any known comet. Trails t1 and t2 represent a much larger integrated brightness than all known cometary trails combined and may therefore be major contributors to the circumsolar dust cloud. We used our Spitzer observations to determine the orbits of these trails and were able to link them to two ("orphan" or type II) trails that were discovered by the Infrared Astronomical Satellite (IRAS) in 1983. The orbits of trails t1 and t2 that we determined by combining the Spitzer and IRAS data have semimajor axes, eccentricities, and inclinations like those of the main-belt asteroids. We therefore propose that trails t1 and t2 were produced by very recent (<~100 kyr old) collisional breakups of small, <~10 km diameter main-belt asteroids

A unified solution for the orbit and light-time effect in the V505 Sgr system

The multiple system V505 Sagittarii is composed of at least three stars: a compact eclipsing pair and a distant component, which orbit is measured directly using speckle interferometry. In order to explain the observed orbit of the third body in V505 Sagittarii and also other observable quantities, namely the minima timings of the eclipsing binary and two different radial velocities in the spectrum, we thoroughly test a fourth-body hypothesis - a perturbation by a dim, yet-unobserved object. We use an N-body numerical integrator to simulate future and past orbital evolution of 3 or 4 components in this system. We construct a suitable chi^2 metric from all available speckle-interferometry, minima-timings and radial-velocity data and we scan a part of a parameter space to get at least some of allowed solutions. In principle, we are able to explain all observable quantities by a presence of a fourth body, but the resulting likelihood of this hypothesis is very low. We also discuss other theoretical explanations of the minima timings variations. Further observations of the minima timings during the next decade or high-resolution spectroscopic data can significantly constrain the model

Yarkovsky-O’Keefe-Radzievskii-Paddack effect on tumbling objects

A semi-analytical model of the Yarkovsky–O’Keefe–Radzievskii–Paddack (YORP) effect on an asteroid spin in a non-principal axis rotation state is developed. The model describes the spin-state evolution in Deprit–Elipe variables, first-order averaged with respect to rotation and Keplerian orbital motion. Assuming zero conductivity, the YORP torque is represented by spherical harmonic series with vectorial coefficients, allowing us to use any degree and order of approximation. Within the quadrupole approximation of the illumination function we find the same first integrals involving rotational momentum, obliquity and dynamical inertia that were obtained by Cicaló & Scheeres. The integrals do not exist when higher degree terms of the illumination function are included, and then the asymptotic states known from Vokrouhlický et al. appear. This resolves an apparent contradiction between earlier results. Averaged equations of motion admit stable and unstable limit cycle solutions that were not previously detected. Non-averaged numerical integration by the Taylor series method for an exemplary shape of 3103 Eger is in good agreement with the semi-analytical theory